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Aberdeen  Proving  Ground,  Md. 
March  I962 


AN  APPLICATION  OF  THE  JET-FORMATION  THEORY 
TO  A  105M4  SHAPED  CHARGE 


ABSTRACT 

The  eq.uations  describing  the  cone-collapse  and  jet-formation  process 
are  shovn  to  be  in  excellent  agreement  with  the  emerging- jet-velocity, 
radioactive-tracer,  and  sli;g- recovery  data  for  the  105n!m  unconfined  test 
charge.  Because  the  number  of  parameters  measured  experimentally  exceeds 
the  number  needed  to  solve  the  equations,  compatibility  of  the  data  with 
the  equations  shows  that  the  non-steady  theory  correctly  describes  the 
formation  of  that  portion  of  the  jet  able  to  penetrate  mild  steel. 

A  polynomial  representation  of  the  collapse  velocity  is  obtained  which 
is  believed  sufficiently  accurate  for  most  problems  requiring  values  for 
this  parameter. 
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INTRODUCTION 


This  report  presents  numerical  results  from  an  accurate 
computation  of  the  liner  collapse  velocity  for  a  105™i  shaped  charge. 

A  drawing  of  the  charge,  together  with  a  table  of  pertinent  dimensions, 
is  presented  in  Figure  1.  In  addition  to  obtaining  values  for  the 
collapse  velocity,  it  is  also  shown  that  the  available  data  for  this 
charge  are  compatible  with  the  currently  accepted  theoiy  of  liner 
collapse  and  jet  formation. 

The  first  published  theory  describing  the  formation  of  a  high-speed 

metallic  jet  by  a  lined-cavity  charge  appeared  in  the  Journal  of  Applied 

Physics  for  June,  1948^.  The  mathematical  description  of  the  process 

was  based  on  steady-state  hydrodynamics  and  predicted  the  formation 

of  a  jet  whose  length  was  constant  and  equal  to  the  slant  height  of  the 

cone.  There  was,  however,  ample  experimental  evidence  to  show  that  jets 

formed  by  shaped  charges  contained  a  spectrum  of  velocities  with  the 

front  of  the  jet  traveling  much  faster  thau  the  rear.  In  order  to 

2 

describe  the  velocity  gradient  in  the  jet,  Pugh  and  his  associates 
developed  a  more  general  theory  of  the  jet-formation  process.  The 
generalized  theory  is  based  on  the  same  steady-state  hydrodynamic 
concepts  as  the  original  theory;  however,  the  collapse  velocities  of 
the  liner  elements  are  no  longer  assumed  to  be  the  same  for  every 
element,  but  are  allowed  to  vary  according  to  the  position  of  the  liner 
element  within  the  charge.  The  pertinent  equations  describing  the 
collapse  of  a  conical  shaped-charge  liner  by  a  plane  detonation  wave 


propagating  parallel  to  the  axis  of  the  charge  are: 

5  =  sin"^(VQ/2U)  (l) 

V.  =  V  csc(p/2)cos(a  +  6  -  p/2)  (2) 

0  o 

dM  /dM  =  cos^(p/2)  (3) 

s 

sin(ci!  +  25)  -  X  sina  [l  -  tan6tan(a  +  ^)J 

tanp  •  - ^ - - - - -  (4) 

cos  (a  +  25;  +  X  sina  pan(a  +  5)  +  tanSj 
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where  a  is  the  half  angle  of  the  conical  liner;  5,  the  angle  between 

the  direction  an  element  of  the  liner  travels  after  being  struck  by 

the  detonation  wave  and  the  normal  to  the  surface  of  the  liner;  V  , 

the  velocity  at  which  the  liner  travels  towards  the  axis ;  V . ,  the 

J 

velocity  of  the  jet  element  formed;  p,  the  angle  between  the  collapsing 

liner  wall  and  the  axis;  M  ,  the  mass  of  the  slug;  M,  the  mass  of  the 

s 

liner;  and  x,  the  initial  axial  coordinate  of  the  liner  element.  Primed 
quantities  are  obtained  by  differentiating  with  respect  to  x.  The 
velocity  with  which  the  detonation  front  sweeps  along  the  liner  surface 
is  designated  by  U.  For  a  plane  detonation  front  and  a  conical  liner 
U  =  Uj^/cosa,  where  is  the  detonation  velocity  of  the  explosive. 


Equations  (1)  to  (4)  represent  four  independent  relations 
involving  the  five  unknown  quantities,  M  ,  p,  V.,  V  ,  and  6.  In 

S  JO 

principle,  the  experimental  determination  of  any  one  of  the  quantities 

is  sufficient  to  provide  a  complete  description  of  the  process  since 

Equations  (l)  to  (4)  can  then  be  solved  for  the  remaining  four  quantities. 

Furthermore,  an  experimental  check  on  the  accuracy  of  the  theory  can  be 

obtained  by  measuring  any  two  of  the  five  quantities.  The  qualitative 
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accuracy  of  the  jet-formation  theory  was  originally  demonstrated  by 

experimentally  determining  P  (x)  and  V.(x).  Experimental  values  for 

J 

P(x)  were  readily  obtained  from  slug- recovery  experiments  using  (5), 

and  V.(x)  was  determined  from  emerging- jet-velocity  and  jet-collection 
J 

measurements.  By  subtracting  the  slug  mass  from  the  corresponding 

liner  mass,  the  cumulative  mass  of  jet  material  M.(x)  was  determined. 

J 

It  was  then  possible  to  determine  the  value  of  x  corresponding  to  a 

given  amount  of  penetration,  x(P),  by  firing  through  targets  of  various 

thickness  and  collecting  the  unused  jet  material.  The  velocity  of  the 

impinging  jet,  V.(P),  was  obtained  as  a  function  of  penetration  depth 
J 

by  measiiring  the  velocity  of  jets  emerging  from  targets  whose  thicknesses 

were  varied.  The  jet  velocity  corresponding  to  the  liner  element  located 

at  position  x  was  then  obtained  by  eliminating  the  penetration  from  the 

two  functions  x(P)  and  V.(P).  The  method  of  analyzing  the  data  that 

J 

was  used  in  the  original  verification  of  the  jet-formation  theoiy  is 


illustrated  by  the  diagram  shown  in  Figure  2.  The  general  agreement 
between  the  3  calculated  directly  from  slug-recovery  data  and  the  3 
calculated  from  (4)  showed  the  theory  to  be  q.mlitatively  correct. 

It  was  recognized  that  the  original  study  of  the  collapse  process 
could  be  improved  in  at  least  two  respects.  First,  the  jet-collection 
experiment  used  in  determining  V.(x)  was  not  as  accurate  as  desired. 

O 

Second,  the  method  of  analyzing  the  data  required  an  experimental  curve 
to  be  differentiated  in  order  to  determine  3(x),  and  these  values  for 
3(x)  were  used  to  determine  V^.  When  (4)  is  finally  used  to  compute  new 
values  of  3  8-s  a  check  on  the  validity  of  the  theory,  it  is  necessary  to 
obtain  the  derivative  of  V^.  Obviously  this  involves  terms  containing 
the  second  derivative  of  the  original  experimental  data.  For  this  reason, 
there  is  some  question  concerning  the  numerical  accuracy  of  this  method 
for  determining  the  compatibility  of  the  experimental  data  with  (l)  to  (4). 
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Eichelberger  ,  in  a  later  study  of  the  problem,  completely 
rejected  the  jet-collection  data.  He  obtained  by  integrating  (4) 
using  experimental  values  of  3  determined  from  new  and  more  precise 
slug-recovery  experiments.  The  values  of  and  3  were  used  to  compute 
the  Jet  velocity  from  which  it  was  possible  to  obtain  the  spatial 
position  of  the  jet  elements  at  any  given  instant  of  time.  These 
results  were  then  compared  with  experimental  jet-velocity  data  in  order 
to  obtain  a  check  on  the  accuracy  of  the  theory.  The  largest  source 
of  uncertainty  in  this  analysis  arises  from  the  difficulty  in  choosing 
the  initial  conditions  for  integrating  (4). 

Uncertainty  in  the  V.(x)  curve  due  to  inaccuracies  in  the  jet- 
collection  experiment  can  be  circumvented  by  a  more  precise  experiment 
in  which  a  radioactive  tracer  is  plated  on  the  inside  surface  of  the 
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liner.  The  first  experiments  with  radioactive  tracers  were  no  more 
precise  than  the  original  jet-collection  experiments.  However,  it  was 
apparent  that  the  radioactive-tracer  technique  coiild  be  used  to  obtain 
reliable  data  for  large  caliber  charges,  such  as  the  105mm  test  charges 
lised  at  the  Ballistic  Research  Laboratories.  During  the  past  several 
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years,  complete  data  have  heen  chtained  for  unconfined  charges 
containing  10 5mm  conical  copper  liners.  This  report  contains  an 
analysis  of  the  data  obtained  and  represents  a  complete  description 
of  the  collapse  process  for  these  charges. 


REDUCTION  OF  EXPERIMEMTAL  RESULTS 


The  penetration  produced  by  a  specific  zonaJ.  element  was 


110 


on  the  inside  sirrface  of  the 


determined  by  plating  a  band  of  Ag 
liner  * ,  Ag  wB-s  chosen  uecausB  ii  is  rcEdily  dstected^  giving  off 
both  strong  gamma  and  high-energy  beta  radiation.  The  physical 
properties  of  Ag^^^  are  similar  to  those  of  copper  and  it  possesses 


o  ^vi4*  1  nr  T rf  Vkto  T  ^ _ T  ■f  •P.a  /  Orrn  ^ 


I pVs^  T  ^*ntM  rvViQ^rr^ 

-i-W  ^UMXi  W  I  tC4«X 


a  large  hole  in  which  the  transportation  of  radioactive  material  from 
one  part  to  another  by  succeeding  portions  of  the  jet  is  sufficiently 


small  so  that  reproducible  data  can  be  obtained.  The  charges  were 
encased  in  a  0.375  in.  thick  plastic  body  and  fired  from  an  8.25  in. 
standoff  into  stacked  6.0  in.  by  6.0  in.  by  1.0  in.  mild-steel  plates. 
The  results  were  sufficiently  reproducible  to  determine  the  location 
of  the  penetration  produced  by  the  tagged  jet  element  to  within 
0.50  inch.  The  experimental  data  were  interpolated  to  obtain  the 
penetration  depths  corresponding  to  values  of  x  between  4.0  cm  and 
8.0  cm  inclusive.  Values  of  P  were  determined  for  every  0.5  cm  change 
in  X  within  this  range  and  checked  for  numerical  accuracy  by  inverse 
interpolation.  The  results  are  presented  in  Table  I. 


The  emerging  jet  velocities  were  determined  by  firing  the  charges 
from  an  8.25  in.  standoff  into  targets  whose  thicknesses  were  varied 
systematically  .  It  should  be  noted  that  these  charges  were  not  encased 
in  the  light-weight  plastic  bodies  used  in  the  radioactive  work. 

However,  it  was  believed  that  the  0.575  in.  plastic  casing  did  not 
seriously  adter  the  jet-formation  process,  and  all  the  data  were 
regarded  as  representative  of  the  charge  Illustrated  in  Figure  1.  The 
emerging  jet  velocities  were  determined  using  an  oscilloscope  and  a 
system  of -grids,  which  measured  the  time  req.ulred  for  the  jet  to 
travel  between  two  known  points  in  space.  This  technique  for  measuring 
the  emerging  jet  velocities  was  chosen  because  it  is  possible  to  obtain 
reliable  data  at  velocities  below  those  for  which  optical  methods  can  be 
used.  The  data  obtained  with  the  electronic  system  were  checked 
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against  similar  data  obtained  with  a  streak  camera  for  those  velocitio's 
where  the  optical  techniques  were  applicable.  The  experimental  data 


w’ere  interpolated  t 


u  u  u  o 


aln  jet  velocities  coriesiXjnding  to  penetration 


depths  that  were  multiples  of  1.0  in.  or  2.5^  cm.  The  results  were 
smoothed  slightly  to  eliminate  fluctuations  in  the  first  differences 
and  checked  for  numerical  accuracy  by  inverse  interpolation.  The 
smoothed  values  of  P  and  V.,  which  were  used  to  represent  the  emerging- 
Jet-velocity  data,  are  presented  in  Table  II. 


The  slug  mass  produced  by  a  liner  section  was  determined  by  firing 
sectioned  105mm  liners  into  water  and  recovering  that  portion  of  the 
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slug  corresponding  to  the  pertinent  liner  section  .  Two  lOpinm  liners 
were  cut  at  a  known  reference  plane  perpendicular  to  the  cone  axis. 

Two  liners  were  used  in  order  to  compensate  for  the  metal  lost  in  the 
machining  process.  The  upper  section  was  weighed,  and  then  the  two 
sections  were  glued  together.  The  charges  were  fired  into  a  10-in. 
gun  barrel  filled  with  water.  Two  recognizable  slugs  were  recovered 
per  shot,  and  the  one  corresponding  to  the  upper  section  of  the  cone 
■was  weighed.  This  was  done  for  a  number  of  reference  planes,  yielding 
slug  mass  as  a  f\inction  of  cone  mass  and,  simultaneously,  cone  mass 
as  a  f-unction  of  position  along  the  cone  axis.  It  should  be  noted 
that  the  light-weight  plastic  bodies  previously  mentioned  were  not 
used  in  the  slug-recovery  program. 


The  experimental  data  were  interpolated  to  obtain  numerical  values 
of  slug  mass  versus  cone  mass  for  equal  Intervals  of  the  cone  mass. 

The  values  were  smoothed  to  eliminate  fluctuations  in  the  first 
differences  and  checked  for  numerical  accuracy  by  inverse  interpolation. 
Results  are  presented  in  Table  III.  The  cone  mass  was  also  obtained  for 
equal  intervals  of  x  in  the  range  between  4.0  cm  and  8.0  cm,  checked  by 
inverse  interpolation,  and  presented  in  Table  lY.  Since  the  entire 
penetration  is  produced  by  jet  elements  originating  between  4.0  cm  and 
8.0  cm  from  the  apex  of  the  liner,  there  is  no  reason  to  define  the 
function  outside  this  interval  unless  penetrations  into  lowei’  strength 
materials  are  to  be  considered. 
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DETERMINA.TION  OF  THE  COLIAPSE  VELOCITY 

The  TiQn-^^eady  theory  of  jet  formation  discussed  in  reference 
(2)  assvimes  a  plane  detonation  wave  throughout  the  collapse  process. 
This  was  a  v^id  approximation  at ,  the  time,  as  the  charges  were  of 
considerable  len^h.  The  experimental  results  analyzed  here  were 
■obtained  from  short  charges;  hence,  the  plane-wave  approximation  is 
no  longer  valid.  The 'theory  has  been  generalized  for  the  short-charge 
geometry  by  assuming  a  spherical  detonation  wave  emerging  from  the 
point  of  initiation.  The  primary  difference  between  the  equations 
for  a  plane  'wave  and  those  for  point  initiation  is  that  U  is  no  longer 
constant,  but  instead  is  given  by  the  relation 


U  =  Uj^cosy 


(5) 


where,  for  eaph  point  on  the  liner,  7  is  the  minimum  angle  between 
the  normal  to  the  spherical  wave  front  and  the  cone  wall  at  the 
intersection  of  the  detonation  wave  and  the  cone  -wall.  One  now 
designates  the  position  of  the  liner  element  as  a  function  of  time  by 
cylindrical  coordinates  r(t)  and  z(t)  and  determines,  as  in  reference  (2), 
tanp  by  evaluating  (uv/Sz)^  when  the  liner  reaches  the  axis  (r  =  O). 
Equation  (4)  becomes: 


tang  = 


_  A 

OfcLllU  X  V  UUfciH. 

O 


A  OcllM.x\u'  ucau-ui. 


Tri  Ar  ^ 

V  '  /  V  I 

o'  o' 


1  +  X  tetna  Qb’  +  (V^/V^)tanAj  -  T'V^sinA 


(6) 


=  Pyi  A/  _  III  Arl  +sr>a.  /i  =  rv  4-  a.  and  is  the  time  at 

L  o'  o  '  '  J  '  "  -  -  V-  / 

which  the  detonation  wave  reaches  a  position  x  on  the  liner. 


fTTV}^  o  »-»o  T  ir«r- *?  n  •?  r*  ^  1  1  ii  c»+ -no  4*  o  rR  "h^r  H  1  n  £0^5:^  TTl  in  T^’l  CnmA 

The  slug-recovery. data  were  differentiated  once  in  order  to  obtain  a 

first  approximation  for  the.  collapse  angle,  p(x).  The  radioactive- 

tracer  and  emerging- Jet-velocity  data  were  used  to  define  an  empirical 

V.(x).  Estimates  of  the  .collapse  velocity  for  the  first  portions  of 
J  , 

the  jet  were  computed  using  (l)  and  (2).  The  values  of  the  collapse 
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velocity  scattered  substantially  due  to  fluctuations  in  the  derivative 
of  the  slug- recovery  data;  however,  a  value  between  0.20  cra/tisec  and 
0.21  cm/usec  was  indicated  for  the  collapse  velocity  corresponding  to 
the  tip  of  the  jet:  i.e.,  the  liner  element  located  at  x  =  4.0  era. 

A  value  was  chosen  for  the  collapse  velocity  of  the  element  located  at 
X  =  4.0  cm,  and  collapse  velocities  were  computed  for  x  >4.0  cm  by 
propagating  the  solution  using  the  first  two  terms  of  a  Taylor  series; 
i.e., 

V^(x  +  Ax)  =  V^(l  +  V^Ax/V^).  (7) 

To  determine  V’/V  ,  (6)  was  solved  for  V’/V  and  numerical  values 

Q  Q^  O  O 

calculated  using  the  approximate  values  of  p  obtained  from  the  slug- 

recovery  data.  At  this  point  new  values  of  p  were  calculated  by 

solving  (l)  and  (2),  using  the  empirical  V.(x)  curve.  New  values  of 

d 

were  calculated  using  the  new  values  of  P  and  the  process  repeated 
until  no  further  change  occurred  in  p.  This  part  of  the  analysis  is 
an  iteration  process  and  is  designated  in  Figure  5  by  the  heavy  arrows. 
The  collapse  angle  calculated  from  the  jet-velocity  curve  using  the 
iteration  process  was  brought  into  satisfactory  agreement  with  the 
slug-recovery  data  by  making  minor  adjustments  in  the  value  of 
corresponding  to  the  tip  of  the  jet. 

A  check  on  the  numerical  accuracy  was  made  by  fitting  the  resultant 
values  of  with  a  fourth-degree  least-squares  polynomial  and  computing 
dM  /dM  and  V..  This  revealed  a  systematic  difference  which  was 

S  J 

attributable  to  use  of  only  two  terms  in  the  Taylor  series.  The  error 
■vra,s  six  parts  in  the  fourth  figure  at  x  =  4.5  cm  and  four  parts  in  the 
third  figure  at  x  =  8.0  cm.  These  differences  were  sufficiently  small 
to  allow  correction  by  a  minor  adjustment  of  the  coefficients  of  the 
first  two  terms  of  the  polynomial.  The  final  results  of  the  analysis 
are  presented  in  Table  V. 
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When  calculating  the  jet  and  sliig  parameters  from  the  polynomial 

for  it  was  noted  that  had  to  have  four  significant  figures  in 

order  to  obtain  a  reasonable  reproduction  of  V  .  in  the  region  4. 0  cm 

J 

X  5*5  cm  because  the,  computed  values  of  V.  were  quite  sensitive 

J 

to  the  derivative  of  V  .  Increased  accuracy  in  the  derivative  of  V 

o  o 

was  achieved  by  computing  the  collapse  velocity  from  the  empirical 
V.(x)  values  and  using  the  sliig- recovery  data  only  as  a  check  on  the 
final  results.  Of  course-  values  of  presented  in  Table  V  are  not 


accurate  to  four  sigiiificant  figures  T3£CS.Li3S 
experimental  data. 


Graphical  comparisons  of  the  results  of  this  analysis  with  the 

experimental  data  are  shown  in  Figure  4  and  Figure  5*  Figure  4 

compares  calculated  from  the  fourth-degree  polynomial  and  the 

experimental  values  determined  at  Carnegie  Institute  of  Technology*. 

The  experimental  values  are  questionable-  in' the  second  figure;  however, 

they  serve  to  demonstrate  that  the  polynomial  for  is  qualitatively 

correct.  Shown  also  in  Figure  4  is  a  comparison  between  experimental 

values  of  V .  and  those  calculated  from  the  polynomial  for  V  .  The 
j  •  o 

latter  comparison  demonstrates  the  accuracy • with  which  the  polynomial 
represents  the  original  data. 


The  jet  mass,  as  well  as  the  jet  velocity,  is  important  in 

determining  the  diameter  of  the  hole  produced  by  the  jet.  It  is 

therefore  necessary  that  values  of  M.  computed  from  the  polynomial  for 

J 

he  in  accord  with  experimental  obBervatlons.  Figure  5  shows  the 

excellent  agreement  between  experimental  values  of  M.,  which  are 

J 

plotted  as  open  circles,  and  the  M=  versus  M  computed  from  the  polynomial 

J 


for  V 


It  should  be  noted  that  M.  is  defined  as  M  -  M  and  there 

J  s 


for< 


includes  approximately  4  g  of  liner  material  dispersed  as  fine  particles 
ahead  of  the  penetrating  jet. 


R.  W.  Watson  and  K.  R.  Eecker,  Private  Communication, 
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CONCLUSIONS 


Although,  a  large  amount  of  shaped-charge  research  has  heen 
completed,  the  collapse  velocity  remains  a  most  difficult  parameter 
to  determine  accurately.  Theoretical  computations  of  have  not  heen 
refined  to  the  point  where  they  can  he  used  without  serious  reservations 
Direct  experimental  observations  have  provided  a  few  values  of  V^,  hut 
the  numerical  accuracy  of  these  measurements  is  questionable  in  the 
second  figure.  The  systematic  analysis  of  slug-recovery,  emerging- 
jet-velocity,  and  radioactive -tracer  data  is  the  only  method  for 
obtaining  accurate  numerical  values  of  V^.  The  data  obtained  with 
105mm  unconfined  test  charges  have  been  carefully  analyzed  and  a 
polynomial  representation  of  the  collapse  velocity  has  been  obtained. 
Because  the  slug- recovery  data  were  used  only  indirectly  in  determining 
the  polynomial  for  V^,  the  excellent  agreement  illustrated  in  Figure  5 
shows  the  non-steady  theory  of  Jet  formation  to  be  adequate  for  that 
portion  of  the  Jet  able  to  penetrate  mild  steel. 
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TABLE  I 


Tabulated  values  of  the  penetration  P  as  a  function  of  the  axial 
position  X  of  the  liner  element  producing  the  penetration.  The 
tabulated  values  are  based  on  radioactive-tracer  data  and  the 
numerical  accuracy  is  checked  by  inverse  interpolation. 
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Inverse  Interpolation 
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TABLE  II 


^  x» 


TabuJLated  vaiijes  of  the  impinging  jet  velocity  V .  as  a  funeuion  oi 

J 

the  penetration  depth  P.  The  tabulated  values  are  based  on  emerging- 
jet-velocity  data  and  the  numerical  accuracy  is  checked  by  inverse 
interpolation i 


Emerging  Jet  Velocity 
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p 

cm 

V  • 

cm/ u-sec 

/j 

cm/ p- sec 

0 

0.701 

0.701 

2.54 

0.643 

0.647 

c:.  /*\Q 

0 . 601 

0.601 

10.16 

0.523 

0.524 

15.24 

0.462 

0.461 

20.52 

o.4o8 

0.409 

27.94 

0.350 

0.549 

35.56 

0.293 

0.295 

43.18 
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18 


TABLE  III 


Tabulated  values  of  the  cumulative  sliog  mass  as  a  function  of 
cumulative  cone  mass.  The  tabulated  values  eire  based  on  slug- 
recovery  data  and  the  numerical  accuracy  is  checked  by  inverse 
interpolation. 
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TABLE  IV 


Tabulated  values  of  ciomulative  cone  mass  as  a  function  of  the  axial 
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values  are  based  on  slug-recovery  data  and  the  numerical 
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checked  by  inverse  interpolation. 
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TABLE  V 


Tabulated  values  of  the  jet  and  slvig  paraineters  computed  from  the 
fourth-degree  polynomial  representation  of  the  collapse  velocity. 
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Figure  2:  The  methcxl  used  to  analyze  the  slug 
recovery  and  jet-collection  data. 
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Figure  5:  Method  used  to  analyze  the  data  obtained 
with  the  lO^mm  unconfined  charge. 
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X-CM 

Figure  4-:  Results  obtained  from  an  analysis  of  the  emerging- jet-velocity  and 
radioactive -tracer  data  for  the  105mm  unconfined  shaped  charge.  The  open 
circles  represent  results  obtained  from  a  numerical  analysis  of  the  data, 
the  triangles  represent  experimental  values  obtained  at  C.  I.  T.,  and  the  curves 
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